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Abstract 

We examine how constraints can be placed on the neutrino component of dark 
matter by an accurate measurement of neutrinoless double beta decay and the 

solar oscillation amplitude. We comment on the alleged evidence for Qv(3(3 decay. 



The detection of neutrinoless double beta decay would imply the violation of lepton 
number conservation. The process could be induced by Majorana neutrino mass terms, or 
by less trivial modifications of the standard model. Here we consider the former possibility 
wherein there are exactly three left-handed neutrino states with Majorana masses 0]. The 
measurement of 0^(3(3 decay, together with what has been learned from studies of solar 
and atmospheric neutrinos, has direct consequences on the spectrum of neutrino masses and 
therefore on the effects of neutrinos on structure formation. We define what will be necessary 
to determine the neutrino component of dark matter from terrestrial experiments. 

The charged- current eigenstates are related to the mass eigenstates by a unitary trans- 
formation 
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where and q are the sines and cosines of 9i, and V is the diagonal matrix (1, e'", e**^ 2 +°>y 
In Eq. (|l|), 02 and 03 are additional phases for Majorana neutrinos that are not measurable 
in neutrino oscillations; if CP is conserved, the phases in UV are either or vr. 

The solar neutrino data favor the Large Mixing Angle solution with 0.6 < sin^ 26^ < 0.98, 
and 2 x 10"^ eV^ < < 4 x 10"'' eV^ at the 3(7 C.L. [0. Atmospheric neutrino data imply 
sin2 20i > 0.85 and 1.1 x lO^^ eV^ < < 5 x lO'^ eV^ at the 99% C.L. ||. The CHOOZ 
reactor experiment imposes the constraint sin^2^^2 ^ 0.1 at the 95% C.L. and A^ 

are the mass-squared differences relevant to solar and atmospheric neutrino oscillations, 
respectively. 

We choose the mass ordering mi < 1112 < with rrii non-negative. There are two 
possible neutrino mass spectra: 
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where in either case A^ ^ A^ in accord with the previously described experimental data. 
For the normal hierarchy (Case I), mixing is given by Eq. (|I]). The limit on 62 implies that 
there is very little mixing of Ve with the heaviest state. In Case I solar neutrinos oscillate 
primarily between the two lighter mass eigenstates. For the inverted hierarchy (Case II), 



solar neutrinos oscillate primarily between the two nearly degenerate heavier states. In this 
case the mixing is described by interchanging the roles of mi and m^. With a mixing matrix 
obtained from Eq. (|l]) by interchange of the first and third columns of UV, the parameters 
governing neutrino oscillations {6i and 6) retain the same import as those in Case I. The 
limit on 62 again implies that for Case II there is very little mixing of z/g with the lightest 
state. 

The rate of OujSjS decay depends on the magnitude of the Pg-i^e element of the neutrino 
mass matrix which is 

Mee = \clclmi + clslm2e"''^ + slmse"'"\, (Case I) , (4) 
= \clclms + clslm2e''''^ + slmie'''"\, (Case II) . (5) 

The masses rrii may be determined from the lightest mass mi and the mass-squared differ- 
ences. Since the solar mass-squared difference is very small it can be ignored; then setting 
mi = m and = A, 



1712 = ^ 1 "^3 = \nn? + A , (Case I) , (6) 
1712 = = Vm^ + A , (Case II) . (7) 

The lightest mass is related to the sum of neutrino masses (S = Sm^) via 



E = 2m + Vm2 + A , (Case I) 
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+ 2Vm2 + A , (Case II). (9) 



For a given value of M^e, the minimum possible value of m is obtained if the three 
contributions to Mee are in phase, i.e., 02 = 03 = 0. Thus 



Mee4 - si J Ml + A cos 202 

mmin = 7^ , (Case I) , (10) 
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The maximum possible value of m is obtained if the the two smaller contributions to Mee 
are out of phase with the largest contribution {i.e., 02 = 03 = tt when C3 > S3). Then 



Mee4\ COS 2^3 I + S^J M^ + A(c| COs2 2^3 - s|) 

"^max = 4 2oa 4 ' (CaSC I) , (12) 
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The allowed ranges for E are determined from Eqs. and (^. Because 62 is small (sin^ 2^2 < 
0.1 or si < 0.026), its value does not significantly affect the result. (We have confirmed this 
result numerically). The limits on S (for 62 = 0) are 



2Mee + ./M2 ±Acos2 26. 



2M. + V^^±A< S < ^-^^^ , (14) 

where the plus sign applies to the normal hierarchy and the minus sign to the inverted 
hierarchy. The bounds depend on only two oscillation parameters: the scale of atmospheric 
neutrino oscillations (A) and the amplitude of solar neutrino oscillations (sin^2^3). 

Figure |l| shows the allowed bands for S and M^e with several possible values of ^3 within 
its 3cr allowed range. We have fixed A = 3 x 10^'^ eV^ and 62 = 0. The solid line is the ^3- 
independent lower bound on S from Eq. ([1^ . The several dashed lines (labelled by sin^203) 
are upper bounds on S. The bands between solid and dashed lines are the (6^3 -dependent) 
allowed domains of E and M^e- The bands in Fig. 0b terminate at Mee = 0.055 eV because 
Mee > VA for Case II. It is possible that CP violation in the neutrino sector is absent or 
negligible. In this case, the point defined by S and Mee niust lie on one of the bounding 
lines of the allowed region. 

If the recent evidence that 0.05 eV < M^e < 0.84 eV at the 95% C.L. [I is borne out, 
this would imply that 0.1 eV ^ S < 20 eV, which using Qi^h"^ = S/(93.8eV) translates to 

0.001 < fi^/i^ < 0.2, (15) 

where is the fraction of the critical density contributed by neutrinos and h is the di- 
mensionless Hubble constant {Hq = lOO/i km s~^ Mpc~^). CMB measurements and galaxy 
cluster surveys already constrain S to be smaller than 4.4 eV [Vtyh"^ ^ 0.05) at the 95% 
C.L. [[7|. Data from the MAP satellite should either determine S or tighten this constraint 
to about 0.5 eV in the near future 0. A more stringent upper bound on S from terres- 
trial experiments must await the precise determination of ^3 (such as is anticipated from 
KamLAND [110]), and a firmer measurement (or constraint) on M^e ||11||. 



Acknowledgements. This work was supported in part by the NSF under grant No. NSF-PHY- 
0099529, in part by the U.S. Department of Energy under grant Nos. DE-FG02-91ER40676, 
DE-FG02-95ER40896 and DE-FG02-01ER41155, and in part by the Wisconsin Alumni 
Research Foundation. 



4 



References 



[1] For related work see e.g., H. V. Klapdor-Kleingrothaus, H. Pas and A. Y. Smirnov, 
Phys. Rev. D 63, 073005 (2001) [|arXiv:hep-ph/0003219|| ; V. Barger and K. Whisnant, 
Phys. Lett. B 456, 194 (1999) ||arXiv:hep-ph/9904281|| ; W. Rodejohann, Nucl. Phys. B 
597, 110 (2001) ||arXiv:hep-ph/0008044|| ; S. Pascoli, S. T. Petcov and L. Wolfenstein, 
Phys. Lett. B 524, 319 (2002) ||arXiv:hep-ph/0110287|| ; H. Pas and T. J. Weiler, Phys. 
Rev. D 63, 113015 (2001) ||arXiv:hep-ph/0101091|| . 



[2] 
[3] 
[4 



[5 



[6, 

[7: 



[9 
[10 



J. N. Bahcall, M. C. Gonzalez- Garcia and C. Pena-Garay, |arXiv:hep-ph /0111150 



T. Toshito [SuperKamiokande Collaboration], |arXiv:hep-ex/0 105023 . 

M. Apollonio et al. [CHOOZ Collaboration], Phys. Lett. B 466, 415 (1999) [|arXiv:hep- 



ex/99070371| . 



M. Doi, T. Kotani, H. Nishiura, K. Okuda and E. Takasugi, Phys. Lett. B 102, 323 
(1981). 

D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15, 1 (2000). 



X. Wang, M. Tegmark and M. Zaldarriaga, |arXiv:astro-ph/0 105091 



H. V. Klapdor-Kleingrothaus, A. Dietz, H. L. Harney and 1. V. Krivosheina, Mod. Phys 
Lett. A 16, 2409 (2002) ||arXiv:hep-ph/0201231|| . For critiques see: C. E. Aalseth et al. 
arXiv: hep-ex/02020 18| ; F. Feruglio, A. Strumia and F. Vissani, |arXiv:hep-ph/020129l" 



D. J. Eisenstein, W. Hu and M. Tegmark, |arXiv:astro-ph/9807130 . 

V. Barger, D. Marfatia and B. P. Wood, Phys. Lett. B 498, 53 (2001) ||arXiv:hep- 



ph/0011251 



[11] A. Alessandrello et al. [CUORE Collaboration], |arXiv:hep-ex/0201038| ; M. Danilov et 
al, Phys. Lett. B 480, 12 (2000) ||arXiv:hep-ex/0002003|| ; H. Ejiri, J. Engel, R. Hazama, 
P. Krastev, N. Kudomi and R. G. Robertson, Phys. Rev. Lett. 85, 2917 (2000) 
arXiv:nucl-ex/9911008|[ ; H. V. Klapdor-Kleingrothaus et al. [GENIUS Collaboration], 



arXiv:hep-ph/9910205 



5 



A = 3x10"^eV2, 02 = 



100 



10 



> 

CD 



0.1 



0.01 



a) Case I 



tritium p decay 

1 




0.001 
100 



0.01 



0.1 



10 



> 



0.1 



0.01 
0.001 



b) Case II 



10 




Figure 1: S versus Mgg for the a) normal hierarchy and b) inverted hierarchy. The sohd 
hne is the 6'3-independent fower bound on S and the dashed hues are the upper bounds for 
different values of 6^. For the inverted hierarchy, Mee > -\/A. The 95% C.L. bounds from 
tritium beta decay and cosmology [0 are shown. 
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